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ABSTRACT 

The  acoustic  volume  reverberation  in  sediments  of  different  grain  sizes  was 
determined  in  the  laboratory.  Three  sediments,  fine  sand,  aquarium  sand  and  aggregate, 
ranging  in  average  grain  size  from  0.3  mm  to  5.3  mm,  were  selected.  An  additional 
material,  glass  beads,  with  an  average  particle  size  of  1.7  mm,  was  also  studied  to 
investigate  the  influence  of  inhomogeneities  in  particle  composition. 

In  the  experiment  ten  sets  of  50  scattered  pulses  of  pulse  duration  88  jis  and  180 
kHz  carrier  frequency  were  averaged,  and  the  envelope  of  the  decaying  tail  was  fitted 
with  a  logarithmic  curve.  The  decay  constants  for  the  various  sediments  showed  a 
variation  with  the  grain  size  expressed  by  a  =  10.3  -  1(1/0.0035  where  a  is  the  grain 
diameter  in  mm  and  r\  the  decay  constant  in  Np/^is.  The  value  of  a  can  be  estimated 
within  ±  2.8  mm  for  a  decay  constant  of  2.5  x  lO'-^  Np/^s.  This  relation  is  limited  to 
the  specific  equipment  used  in  the  experiment  and  to  the  frequency  of  180  kHz. 
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I.  INTRODUCTION 

The  first  use  of  sound  to  detect  underwater  objects  by  echo-ranging  was  in  World 
War  I,  but  it  was  not  until  World  War  II  that  sonar  development  became  significant. 
Detection  and  ranging  are  the  main  functions  of  sonar  systems,  which  include  the 
fathometer,  a  widely  used  and  very  important  application.  The  fathometer  measures 
the  round-trip  travel  time  of  short  sound  pulses  between  a  transmitter  located  at  the 
water's  surface  and  the  bottom.  If  the  sound  speed  in  the  water  is  known,  depth 
measurements  are  possible  and  the  safety  of  navigation  in  near-shore  areas  is 
improved. 

The  fundamental  principle  is  simple  but  the  process  involves  some  complexity 
due  to  the  variation  in  sound  speed  and  attenuation  in  the  water.  These  processes  are 
inherent  to  the  water,  but  when  a  target  or  the  bottom  is  present,  reflection,  absorption 
and  incoherent  scattering  at  the  interface  between  the  target  or  bottom  must  be 
considered. 

This  is  the  situation  present  in  a  hydrographic  survey  operation:  A  transducer  in 
the  water  sends  sound  pulses  to  the  bottom  which  are  reflected  back.  Along  the  sound 
path  transmission  losses  due  to  spreading,  absorption  and  scattering  occur.  Changes  in 
temperature,  salinity  and  pressure  along  the  path  of  a  sound  ray  cause  refraction, 
which  is  not  a  very  important  problem  in  hydrography,  because  the  transmitter  points 
downward.  Upon  reaching  the  bottom,  sound  is  both  reflected  and  transmitted  into 
the  bottom  so  that  its  amplitude  is  decreased.  Some  of  the  sound  transmitted  into  the 
bottom  is  scattered  backward  where  it  combines  with  the  reflected  pulse,  distorting  the 
initial  pulse  shape.   This  pulse  distortion  is  the  subject  of  this  paper. 

In  hydrography  it  is  the  front  edge  of  the  received  pulse  which  is  normally  of 
interest,  since  it  gives  the  distance  to  the  bottom.  In  addition,  however,  more 
information  may  be  obtained  from  the  same  received  pulse.  In  fact,  the  tail  of  the 
received  signal  can  give  information  about  the  nature  of  the  bottom.  There  is  thus  the 
possibility  to  remotely  classify  the  bottom  type,  an  operation  essential  to  the 
preparation  of  nautical  charts. 

The  intent  of  this  project  is  to  verify  a  laboratory  technique  previously  studied 
and  estabUshed  by  Diaz  (1986)  and  Chang  (1986),  in  which  the  trailing  edge  of  the 
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renected  pulse  is  characterized  in  terms  of  a  straight  line  fitted  to  the  envelope  by  a 
least  square  method.  In  this  study  the  same  experimental  procedures  are  applied  to 
natural  sediments  and  artificial  matenals  to  test  the  theory  relatmg  the  type  of  bottom 
sediment  to  the  received  echo  level. 
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II.  THEORETICAL  CONSIDERATIONS 

A.      DEFINITION  OF  REVERBERATION 

In  any  echo-ranging  or  listening  device  the  desired  signal  is  always  received  in  the 
midst  of  a  certain  amount  of  extraneous  noise  (National  Defense  Research 
Commitee,1969).  Particularly  in  underwater  acoustics,  this  noise  can  arise  from 
different  sources  such  as  breaking  waves,  snapping  shrimp,  surf  and  shipping,  which 
combine  to  produce  broadband  ambient  noise  (Kinsler,  Frey,  Coppens  and  Sanders, 
1983). 

However  if  dealing  with  an  active  system  like  the  sounding  equipment  used  in 
hydrography,  another  type  of  interference  arises,  which  is  directly  due  to  the 
transmission  of  the  pulses  into  the  water.  This  component  is  called  reverberation  and 
results  from  the  existence  of  particles  or  discontinuities  in  the  path  of  the  pulse  that  act 
as  small  scattering  centers,  reradiating  the  acoustic  energy  and  returning  a  weak  echo 
to  the  transducer.  The  combination  of  all  these  small  echoes  makes  the  reverberation. 
Like  ambient  noise,  reverberation  is  usually  undesired,  because  it  masks  the  reflected 
signal. 

1.  Surface,  Bottom  and  Volume  Reverberation 

When  small  scatterers  or  discontinuities  are  observed  along  a  reflecting 
surface,  like  the  water-air  interface,  surface  reverberation  occurs.  It  is  important  when 
a  reflecting  surface  has  a  certain  roughness  compared  to  the  wavelength  of  the  sound. 
When  the  reflection  is  produced  at  the  bottom  surface,  and  if  this  bottom  surface 
shows  some  irregularities  or  small  objects,  then  bottom  reverberation  is  produced.  But, 
if  the  reverberation  is  due  to  scatterers  existing  in  the  body  of  the  propagation  medium, 
i.e.  the  water  column,  then  we  have  volume  reverberation.  In  this  project  the  volume 
reverberation  produced  in  the  sediment  is  of  interest  because,  as  it  was  pointed  out 
before,  the  amount  of  volume  reverberation  observed  can  give  an  indication  of  the  type 
of  sediment  present.  This  suggests  the  possibility  of  rapidly  determining  the  types  of 
bottom  sediments  present  during  hydrographic  soundings. 

2.  Reverberation  Level 

If  an  active  system  with  an  intensity  I^  at  the  transmitter  is  considered,  with  a 
target  at  range  r,  the  intensity  level,  IL(r),  of  the  signal  reaching  the  target  is 
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IL(r)  =  10  log  [I(r)/Ij.gf]  dB  (2.1) 


where 


I(r)  =  Iq/i^  (2.2) 


for  transmission  loss  by  spherical  spreading. 

After  being  scattered  by  a  volume,  dV,  the  intensity  level  depends  on  the  ability  of  the 

target  to  reflect  the  incident  energy.    This  ability  is  expressed  by  S^    the  scattering 

strength. 


S^  =  10  log  s^  (2.3) 


where  s    is  the  scattering  coefficient  per  unit  volume 


Sy  =  <slAn  (2.4) 

and  <y  is  the  backscattering  cross  section  of  the  unit  volume.    Thus,  the  scattered 
intensity,  \^{r)  at  1  m  from  the  target  is 


13(1)  =  (Io;r2)s^V  (2.5) 


where  V  is  the  volume  occupied  by  scatterers  that  contribute  to  the  reverberation 
received  at  the  receiver  at  time  t.  One  essential  step  in  obtaining  the  reverberation 
level  (RL)  is  to  compute  the  scattering  volume  which  is  a  function  of  the  transmitting 
pulse  length  (t),  and  the  directivity  of  the  transducer  (here  assumed  monostatic).  The 
geometry  occurring  in  hydrographic  surveying  is  shown  in  Figure  2.1. 
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Figure  2.1     Echo  sounding  geometry  and  elementary  portion 
of  volume  reverberation  contribution. 
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The  thickness  of  the  reverberating  volume  must  be  such  that  all  the  scatterers  in  it 
contribute  reverberation  at  the  same  time  at  the  receiver,  (Kinsler,  Frey,  Coppens  and 
Sanders,  1983)  and  that  depends  directly  on  the  pulse  length  or 


L  =  ct;2  (2.6) 


where  t  is  the  pulse  duration.    For  a  geometric  cross  section  of  the  beam  pattern,  A, 
we  have 


V  =  Ac(T/2)  (2.7) 


Substituting  this  into  Equation  2.5  gives 


13(1)  =  (lQ/r2)s^c(t/2)  (2.8) 


and  the  intensity  at  the  receiver  becomes 

I3  =  (Io;r'^)s^c(t/2)  (2.9) 

The  volume  reverberation  level  in  the  form  of  a  sonar  equation  is 

RL  =  SL-  40  log  r  +  S^  +  10  log  V  (2.10) 


Urick  (1983)  suggests  that  in  reverberation  backgrounds  it  is  easy  to  fmd  the  echo-to- 
reverberation  ratio  for  a  target  of  known  target  strength: 


EL  =  SL-  2TL  +  TS  (2.11) 
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where 

EL  =  echo  level 

SL  =  source  level 

TS  =  target  strength 

Subtracting  Equation  2.10  from  Equation  2.11  we  obtain 

EL- RL  =  TS -(Sy  +  lOlog  V)  (2.12) 


but  (S^.  +  10  log  V)  is  the  reverberation  target  strength,  and  the  echo-to-reverberation 
ratio  is  only  the  difierence  of  two  strengths,  one  for  the  reflected  echo  and  the  other  for 
the  scattered  one. 

How  can  this  suggestion  be  applied  to  our  project  ?  In  it  we  are  interested  in 
finding  the  RLs  produced  by  different  materials.  Considering  the  reflection  only  at  the 
fluid- sediment  interface  (no  incoherent  surface  or  volume  scattering),  the  reflected 
signal  will  have  the  same  shape,  no  matter  which  material  is  below  that  interface,  being 
reduced  in  amplitude  only  by  the  reflection  coefficient  of  the  sediment. 

If  the  incoherent  scattering  produced  by  the  sediment  particles  is  added  the 
received  scattered  signals  will  differ  in  shape,  (mainly  in  the  trailing  edges)  for  different 
sediments.  These  differences  are  due  to  the  incoherent  scattering  represented  in  the 
sonar  equation  by  RL. 

There  is  no  way  to  measure  directly  the  pure  (coherent)  pulse  reflected  at  the 
water- sediment  interface  because  the  reflected  signal  is  "contaminated"  by  incoherent 
scattering,  but  it  can  be  approximated  by  the  pulse  reflected  from  the  water-air 
interface.  The  only  difference  between  these  is  the  voltage  amplitude  because  the 
specific  acoustic  impedances  are  different,  but  this  can  be  corrected  by  normalizing  the 
received  voltages. 

In  Chang  (1986)  this  pure,  reflected  signal  was  obtained  using  a  transducer 
lying  on  the  bottom  and  transmitting  upward  to  the  water-air  interface.  The  signal 
received  is  then  purely  reflected,  assuming  that  no  surface  reverberation  is  produced  at 
the  interface  or  that  its  magnitude  is  negligible  compared  to  the  volume  reverberation 
produced  in  the  sediment.  This  assumption  also  applies  to  the  water-sediment 
interface,  which  must  be  perfectly  level  and  smooth. 
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By  using  the  same  transmitting  conditions,  but  with  the  transducer  located  at 
the  water-air  interface  and  transmitting  downward  to  the  bottom,  a  different  signal  is 
received  back,  in  the  transducer;  this  represents  the  volume  reverberation. 

Geometrically  this  procedure  can  be  illustrated  with  the  pictures  of  the  two 
distinct  waveform  envelopes.  Besides  the  difference  in  amplitude,  the  most  noticeable 
feature  is  that  the  received  waveform  shows  a  slower  decaying  of  its  tail,  where 
reverberation  is  geometrically  expressed.  Mathematically,  the  decaying  tail  can  be 
represented  by 

V^Xt)=  Voe-«^  (2.13) 


for  the  water-air  interface  and 


V^lt)  =  Voe-nt  (2.14) 


for  the  water- sediment  interface,  where  V^  is  normalized  to  the  same  value  in  each  case 
(Chang,  1986).  The  difference  between  these  two  expressions  gives  a  measure  of  the 
amount  of  volume  reverberation  coming  from  that  particular  sediment. 


AV(t)  =  V_(e-^^  -  e'**^)  (2.15) 
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Figure  2.2    Geometric  interpretation  of  volume  reverberation. 
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III.  BACKGROUND 

Many  research  projects  have  been  devoted  to  the  study  of  reverberation  in  the 
ocean.  However,  most  have  been  concerned  with  scattering  produced  in  the  water 
column  and  only  a  few  with  volume  reverberation  in  marine  sediments. 

A.  LITERATURE  REVIEW 

Urick  (1953)  studied  backscattering  from  the  bottom  as  a  function  of  pulse 
length,  frequency  and  grazing  angle  at  several  locations  in  a  harbor  using  a  tiltable 
transducer  mounted  on  a  barge.  The  results  were  expressed  in  terms  of  scattering 
strength  per  unit  area  of  bottom.  Mackenzie  (1961)  made  bottom-reverberation 
measurements  using  two  low  frequencies  and  different  pulse  lengths  as  well  as  different 
grazing  angles.  He  analysed  the  data  to  obtain  scattering  coefficient  for  the  bottom. 
McKinney  and  Anderson  (1963)  made  the  same  measurements  but  for  three  different 
types  of  sediment:  mud,  sand  and  gravel.  They  concluded  that  the  bottom  sediment  is 
a  major  factor  in  the  backscattering  of  sound  and  that  the  knowledge  of  particle  size 
distribution  of  the  sediment  is  useful  in  estimating  the  reverberation  level.  No  model 
was  formulated. 

Recent  relevant  work  (Dodds  [1982,1984],  Breslau  [1967]  and  Meng  [1982])  deals 
with  the  problem  of  using  backscattering  measurements  to  remotely  classify  bottom 
types.  Dodds  used  a  broadband  acoustic  source  (1-10  Khz)  directed  vertically  to  the 
sea  floor  and  obtains  sonograms  from  which  volume  scattering  and  surface  roughness 
parameters  are  obtained.  Meng  (1982)  discussed  the  classification  of  seabed  sediments, 
based  on  the  echo  envelope  of  normally  reflected  high-frequency  sound  pulses.  This  is 
apparently  similar  to  the  present  work.  It  would  be  interesting  to  investigate  his 
procedure,  processing  and  achievements.  However  this  paper  was  published  in 
Chinese,  and  only  a  short  abstract  is  available  in  English. 

B.  VOLUME  REVERBERATION  MODELS 

The  backscattering  of  sound  from  the  ocean  bottom  has  been  studied  for  many 
years,  but  apparently  only  Morse  and  Ingard  (1968)  and  Urick  (1983)  have  developed 
theoretical  models  appliable  to  the  present  study. 
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These  models  have,  however,  a  restricted  application  since  they  make  several 
assumptions  that  are  not  satisfied  in  the  present  case,  where  the  wavelength  is  8  mm 
and  the  average  diameter  of  the  scatterers  is  between  0.3  and  5.3  mm. 

Morse  and  Ingard's  assumptions  and  their  relevance  to  the  present  study  are: 

1)  Scattering  objects  are  small  compared  with  the  wavelength, 
(satisfied  here  only  for  the  fine  sand). 

2)  Dimensions  of  the  scattering  region  are  much  larger  than  the 
wavelength. 

(Assumption  satisfied  for  all  four  sediments). 

3)  The  scatterers  are  all  spheres.  (Not  too  bad  an  assumption  for  all 
four  sediments). 

4)  The  scatterers  are  sparsely  packed  and  no  multiple  scattering  is 
considered.   (Ver\'  poor  assumption  for  all  sediments). 

For  monostatic  sonars  Morse  and  Ingard  give  the  volume  scattering  coefficient 
per  unit  volume  as 


s,.  =  N(2rt)l'Va8(Y^-Yp)^  (3.1) 

where  "a"  is  the  radius  of  the  scatterers,  N  the  number  of  scatterers  per  unit  volume 
and  k  is  the  wave  number, 

'''k  ^  ^'^n"'^o^' *^o 

where  K^^  is  the  compressibility  of  the  scatterer,  Kq  the  compressibility  of  the  fluid  and 

yp  =  (^Pn-^PO)  '  (2Pn-Po) 

where  p^^  is  the  scatterer  density  and  p^  the  fluid  density. 

The  scattering  coefficient  per  unit  volume  as  given  by  Urick  (1983)  is 


s^=N(T  (3.2) 
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where  <T,  the  backscattering  cross  section,  is  given  by 

<T  =  2.8  (71  a^Xk  a)^  (3.3) 

where  "a"  is  the  scatterer  dimension.    This  model  is  good  for  a  single  sphere  (no 
multiple  scattering)  and  for  values  of  Ka  <  <  1. 

C.       DERIVED  MODEL 

The  present  project  is  an  extension  of  the  recent  work  developed  by  Diaz  (1986) 
and  Chang  (1986).  Chang  measured  volume  reverberation  on  aggregate  a.nd  fine  sand 
to  compute  the  corresponding  volume-scattering  coefficients.  Applying  Morse  and 
Ingards  (1968)  and  Urick's  (1983)  models  Chang  did  not  get  satisfactor\'  results, 
verifying  that  these  models  are  not  appropriate  for  ocean  sediments.  A  simple  model 
was  then  set  up: 

The  envelope  of  the  electrical  pulse  shape  can  be  described  by  a  square  wave  of 
pulse  duration  T 


I(t)  =  I  JH(t)  -  H(t-T)]  (3.4) 


where 

I(t)  =  intensity  at  time  t  after  the  beginning  of  the  pulse 

Iq  =  constant  intensity  amplitude 

T  =  pulse  length 

H(t)  =  unit  step  function  shown  in  Figure  3.1  a) 

0  t<0 
H(t)  = 

1  t>0 

H(t-T)  is  then:  (Figure  3.1  b)) 

0  t-T<0 

H(t-T)  = 

1  t-T>0 
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The  transmitted  pulse  is  not  exactly  square  due  to  damping  of  the  transducer,  but  it 
is  tranformed  into  an  acoustic  pulse  with  an  exponential  rise,  a  flattened  top  and  an 
exponential  decay.   This  shape  can  be  described  by:     Figure  3.1  c 


ill)  =  (l-e-«M  H(t)  +  (l-e-*^^)H(t-r) 


(3.5) 


Figure  3.1     Step  unit  function  and  transmitted  acoustic  pulse. 

Where  i  is  the  time  measured  from  the  beginning  of  the  pulse  and  a  the  rise  and  decay 
constants.  For  a  transducer  located  in  the  water  sending  pulses  towards  a  smoothed 
surface  the  intensity  received  by  the  same  transducer  must  be  according  to  Chang 
(1986) 


I  =  [J  Tp  s^  e-n^  Iq  n:t-2^/c)d^'(r  +  ^)2  +  R-I^  /r^l/r^ 


(3.6) 


where: 
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T-  =  intensity  transmission  coefficient 

P  =  absortion  cofTicient  in  the  sediment  =  0.25  f  (Urick,1979) 
^     =  affected  sediment  depth  =  c  T/2 
R-  =  intensity  reflection  coefficient 
r        =  distance  from  transducer  to  the  sediment  surface 
The  first  term  inside  the  square  brackets  represents  the  scattered  component  of  the 
scattered  intensity  for  the  volume  of  the  sand  and  the  second  the  reflected  one. 

For  r  >  >  ^ 


I5  =  (T-2  s^  I^^r^  )  Je-P^  f{i-2^lc)d^  (3.7) 

where 

f]:t-2^/c)  =  [l-e-**(^-2^'^)]H(t-24/c)  +  [l-e-«^^-'^-2^'^)]H(t-T-2^/c)  (3.8) 

After  integration  between  ^  =  0  and  S,  =  cT/2 

Is  =  A  [(l-e-^P^'2)/p-(e-cPt/2.g.at)/(.p  +  2a/c)  (3.9) 

-f  (e-^P(^-^)/2.i)/p  +  (e-^P(^-^)/2.e-«(^-^))/(-P  +  2a/c)] 
where 

A  =  Ti^s^I^/r"^  (3.10) 

With  some  algebraic  operations 
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I3  =  A{l/p[e-*2(,cPT/2.i)j  ^3_jj^ 

and  the  reflected  component 

I,=  IoR-rV'-l)e-«^  (3.12) 


Adding  the  two  components  incoherently,  the  total  (reflected  and 
scattered)  intensity  received  by  the  transducer  will  be 


I  =  B(T-2s^/P[e-'=P^-(ecP^''2.i>,  +  -j-.23^,(.p  +  2a/c)  (3.13) 


[e-^P^Vf^^''--I)  +  e-«Vl-e-«^)]  +  R.(e-«^-l)e-«^} 


where  B  =  I^/r 

The  expression  looks  quite  complex  but  all  the  variables  are  known  except  s^.    If 
the  part  included  inside  braces  is  called  C,  then 

I  =  do  r-^^  C 

The  quantity  inside  brackets  is  the  received  intensity  considering  only  spherical 
spreading  and  C  accounts  for  the  scattered  and  reflected  intensity.  The  only  interesting 
part  is  for  t  >  T,  that  corresponds  to  the  decay  portion.  That  is  also  what  is  measured 
in  the  experiment,  i.e.   the  decay  of  the  pulse  tail. 

To  find  the  value  of  s^  Chang  combined  the  expression  for  C  with  his 
experimental  results.  Guessing  values  for  Sy  the  corresponding  values  for  C  were 
obtained.  He  repeated  then  the  arithmetic  procedure  until  the  value  of  the  log  C  ' 
that  expresses  the  slope  decay  became  close  to  the  decay  constant  obtained 
experimentally.  The  value  guessed  for  s^  is  then  the  volume  scattering  coefficient  for 
the  sediment  studied. 
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This  simple  model  appears  to  approximate  the  volume  scattering.  Further  studv 
of  volume  reverberation  from  ocean  sed.ments  to  improve  the  model  v^as  thus 
suggested. 
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IV.  EXPERIMENT 

A.       LNSTRUMENTATION 

1.  Experiment  Environment 

All  of  the  preliminary  and  main  experiments  were  conducted  in  two  steel- 
bound  glass  tanks  filled  with  firesh  water  and  having  interior  measurements  of  0.7  m  by 
0.7  m  by  0.6  m.    Four  different  materials,  y7/7g  sand,  aquarium  sand,  aggregate 
and  glass  beads,  were  placed  at  the  bottom  of  the  tank  to  serve  as  "sediments". 

The  tank  had  0.35  m  of  water  overlying  a  0.25  m  thick  layer  of  "sediment". 
Specifically,  the  aggregate  and  the  glass  beads  were  kept  inside  a  plastic  bucket  during 
the  experiment  to  reduce  the  amount  of  material  used.  Because  the  transducer  beam  is 
narrow  (projected  spot  is  a  circle  with  0.1  m  radius  at  a  0.3-m  range),  the  beam  fell 
completely  inside  the  bucket  (Figure  4.1). 

The  fine  sand,  also  studied  by  Bradshaw  (1981),  Diaz  (1986)  and  Chang 
(1986),  had  an  average  size  of  0.3  mm  and  varied  from  0.15  mm  to  0.7  mm.  The 
aggregate,  also  used  by  Diaz  (1986)  and  Chang  (1986),  was  less  homogeneous  in  size 
and  composition  than  xhtfine  sand  and  had  an  average  size  of  5.3  mm. 

It  was  convenient  to  make  measurements  using  a  sediment  size  between  0.3 
mm  and  5.3  mm  for  a  better  evaluation  of  the  volume  reverberation  behavior  with  the 
grain  sediment  size,  and  for  that  a  third  sediment  was  selected,  the  aquarium  sand.  This 
one  was  sieved  and  sized  to  an  average  size  of  1.72  mm  and  varied  from  0.84  mm  to 
3.36  mm  with  many  inhomogeneities  in  composition  and  shape.  It  was  expected  that 
these  variations  in  size,  shape  and  constitution  for  the  same  sediment  would  raise  some 
difficulties  and  this  was  confirmed  during  the  experiment.  The  idea  of  performing  the 
experiment  on  an  ideal  homogeneous  sediment  lead  to  the  use  of  spherical  glass  beads, 
whose  density  (2.5  g/cm  )  is  not  far  from  the  real  sediments'  densities  and  size  is  1.7  ± 
0.3  mm. 

The  four  different  materials,  namely  i]\q  fine  sand,  the  aggregate,  the  aquarium 
sands  and  the  glass  beads,  were  kept  in  water  for  two  to  three  days  before 
measurements  were  made.  A  water  jet  from  a  hydraulic  pump,  or  manually  stirring 
was  used  to  remove  trapped  air  bubbles  before  and  during  the  measurement  process. 
Bleach  was  added  to  the  water  to  control  biological  growth,  and  a  wetting  agent  was 
also  mixed  with  water  to  help  remove  air  bubbles. 
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Figure  4.1     nxperimcntal  Geometry  Design. 
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2.  Electronic  Equipment 

A  General  Radio  Model  1310  Oscillator  generated  a  181-kHz  signal  that  was 
simultaneously  fed  into  a  frequency  counter  and  a  General  Radio  type  1396-A  Tone 
Burst  Generator  to  produce  a  16-cycle  pulse,  which  was  then  amplified  by  a  Hewlett 
Packard  467A  power  amplifier  before  passing  to  a  Datasonics  Transmit,' Receive  (T/R) 
switch  that  controlled  the  transmission  and  reception  by  the  same  transducer.  The 
transmitted  signal  was  used  to  trigger  an  analog  oscilloscope  (Tektronix  Type  RM  503) 
and  a  Nicolet  model  3091  Digital  oscilloscope. 

The  transducer,  a  Naval  Research  Laboratory  F-4I  circular  piston  type  with 
an  8.8-cm  diameter  active  face,  has  a  resonance  frequency  of  182  kHz  and  a  half-beam 
width  of  10  degrees  (Diaz,  1986). 

The  signal  received  by  the  transducer  was  amplified  20  dB  by  a  465A  Hewlett- 
Packard  pre-amplifier  and  then  filtered  by  a  Spencer  Kennedy  Laboratories,  Inc.  Model 
302  electronic  filter,  used  as  a  high  pass  filter  to  avoid  lov/  frequency  noises.  The 
filtered  signal  was  finally  passed  to  the  oscilloscope  which  was  controlled  by  a  3421 A 
Hewlett  Packard  Data  Acquisition/Control  unit  attached  to  an  Hewlett  Packard  86 
desk-top  computer. 

All  the  equipment  worked  reasonably  well  :  The  frequency  generator  showed  a 
drift  of  ±  150  Hz  during  a  complete  set  of  10  measurements,  but  this  did  cause  any 
significant  changes  in  the  pulse  waveform,  as  long  as  the  transducer  was  operated  near 
its  resonance.  The  digital  oscilloscope  was  controlled  by  the  computer  and  the  the 
data  acquisition/control  unit,  but  sometimes  during  the  measuring  process  the 
oscilloscope  stayed  locked  in  the  STORE  mode  making  it  necessary  to  interrupt  the 
measurements  and  reset  it. 

Each  set  of  measurements  was  the  result  of  50  samples  and  took  about  45 
minutes.  A  faster  sampling  rate  could  speed  up  the  measuring  process  and  probably 
give  us  a  more  accurate  waveform  reconstitution.  A  block  diagram  of  the  equipment  is 
shown  in  Figure  4.1. 

B.       PRELIMINARY  MEASUREMENTS 

Preliminary  measurements  of  density,  speed  of  sound  and  reflection  coefficient 
were  made  for  later  computations  of  acoustic  characteristics. 


29 


Frequency 
Generator 


Analog 
Oscilloscope 


-) — 


Frequency 
Counter 


Tone  Burst 
Generator 


■"X. 


Power 
Amplifier 


T/R 
Switch 


High  Pass 
Filter 


Digital 
Oscilloscope 


Computer 


F-41  Transducer 


Pre-amplifier 


TrrTTT-rrrm-m 


Figure  4.2     Hlcctronic  Equipment  block  diagram. 


30 


1.  Density 

Bradshaw  (1981)  measured  the  water-saturated  density,  p^^^'  ^^^  ^^  density,  p,  of  fine 

sand  and  obtained: 

p^^  =  1.98  ±  0.03  g/cm^ 

p  =  2.69  ±  0.01  g;cm^ 

For  the  aggregate  Diaz  (1986)  obtained  : 

D    ■     =  1.97  2'cm 
*^rmx       *       s^>-iii 

p  =  2.64  g/cm-' 

Density  measurements  were  also  performed  for  the  aquariumsand  and  for  the  glass 

beads.   The  results  are  summarized  in  Table  1. 

2.  Speed  of  Sound 

The  measurement  of  sound  speed  in  the  sediment  was  performed  inside  the 
glass  tank  using  three  LCIO  transducers,  one  acting  as  a  transmitter  and  the  other  two 
as  receivers.  The  transducers  were  buried  in  the  sediment  at  the  same  depth.  The 
transmitter  and  the  closest  receiver  were  kept  fixed  and  the  third  transducer  was  moved 
to  different  positions  along  the  line  connecting  the  first  two  transducers. 

The  difTerence  in  reception  time  at  the  two  receivers  was  measured  by  means 
of  the  digital  oscilloscope  (Figure  4.3),  and  by  measuring  the  distance  between  the 
receivers  it  was  possible  to  calculate  the  sound  speed.  Fifteen  measurements  were 
made  on  the  aquarium  sand  resulting  in  c  =   1588  ±  22  m/s. 

Hamilton  et  al.  (1956)  studied  the  physical  properties  of  marine  sediments 
from  off  the  coast  of  San  Diego.  California,  and  established  some  relations  between  the 
sediment  grain  sizes  and  their  physical  properties  such  as  porosity,  density,  speed  of 
sound,  etc.  They  observed  a  general  increase  of  sound  speed  with  increase  in  medium 
grain  size.  Thus,  the  value  of  the  sound  speed  for  the  aquarium  sand  was  smaller  than 
the  1610  m/s  obtained  by  Bradshaw  (1981)  for  Iht  fine  sand,  and  does  not  agree  with 
Hamilton  et  al..  the  last  conclusion.  To  investigate  the  consistency  of  that  value,  a 
second  set  of  measurements  was  made  and  the  result  was  c  =  1590  ±  42  m;s, 
confirming  the  first  one. 

The  same  experimental  procedure  was  conducted  on  the  glass  beads  and  the 
results  were  c  =  1822  ±  48  m/s 
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Figure  4.3     Detail  of  speed  of  sound  measurement  procedure. 

3.  Reflection  Coefficient 

The  reflection  coetTicient  for  sound  rays  passing  from  one  medium  to  another 
is  the  ratio  between  the  rcOectcd  and  the  incident  acoustic  pressures,  and  for  normal 
incidence  it  depends  on  the  specific  acoustic  impedance  pc,  of  each  medium.  The 
reflection  coefiicients  for  aquarium  sand  and  glass  beads  were  measured  as  the  ratio 
between  the  maximum  voltage  amplitude  in  the  signal  reflected  from  the  sediment  and 
the  signal  reflected  by  the  water-air  interface  at  the  same  distance. 


R  =  v^/Vi 


(4.1) 


The  reflection  coefTicienl  can  also  be  computed  using 


R  =  (r2-  rj)'(r2  +  rj) 


(4.2) 
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where 

rj  =  pjCj   is  the  specific  acoustic  impedance  for  the  water 

r-)  =  92^2    ^^^  ^^^  sediment 

The  measured,  R^^^,  and  computed,  R^,  values  of  the  reflection  coefficient,  as 
well  as  other  physical  characteristics  are  summarized  in  Table  1. 


TABLE  1 

SEDIMENT  PHYSICAL  CHARACTERISTICS 

fine                aggregate 
sand 

aquarium 
sand 

glass 
beads 

a 

0.3                      5.3 

1.72 

1.7 

(mm) 

(±  0.3)              (±  0.51) 

(±  0.24) 

(±  0.3) 

P 

2.69                     2.64 

2.59 

2.50 

(g/cm^) 

(±  0.01)            (±  0.05) 

(±  0.05) 

{ ±  0.05) 

Pmix 

1.98                     1.97 

2.02 

1.91 

(g;cm^) 

(  ±  0.03)            {  ±  0.05) 

(±  0.02) 

(±  0.01) 

c 

1610                    1555 

1588 

1822 

(m,'s) 

(±  22) 

(±  48) 

Rm 

0.36                    0.36 

0.36 

0.53 

(±  0.09) 

(±  0.05) 

Re 

0.38                    0.52 

0.38 

0.41 

(±  0.01) 

(±  0.01) 
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C.       MEASUREMENT  PROCEDURES 

1 .  Experiment  Technique 

To  measure  the  exponential  decay  of  the  tails  of  the  pulses  scattered  from  the 
sediment  the  experimental  procedure  of  Diaz  (1986)  and  Chang  (1986)  was  used.  The 
transducer  was  clamped  to  a  horizontal  wooden  rail  that  kept  the  transducer  level  with 
its  active  face  pointing  downward  to  the  sediment  as  shown  in  Figure  4.1. 

In  his  experiment  Chang  observed  large  variations  in  the  received  pulse  shapes 
between  different  samples  of  the  same  sediment  and  concluded  that  the  variations  were 
due  to  inhomogeneities  in  the  sediment.  Because  of  this,  he  did  not  obtain  consistent 
results.  To  smooth  the  variations  in  the  results  a  set  of  ten  measurements  was 
performed  with  various  transducer  positions,  and  the  results  were  averaged 

For  the  aggregate  sand  it  is  not  difficult  to  accept  that  important 
inhomogeneities  are  present  in  its  components,  but  the  fine  sand,  looks  quite 
homogeneous.  However,  the  waveforms  did  show  large  variations.  By  performing  the 
experiment  using  a  thin  metal  plate  covering  the  sediment,  Chang  showed  that  those 
inhomogeneities  were  interior  to  the  sediment. 

To  investigate  these  facts  further,  the  first  experiment  performed  in  this 
present  work  consisted  of  keeping  the  transducer  in  a  fixed  position  and,  between  each 
measurement  the  sediment  was  stirred  using  a  water  jet  produced  by  a  pump.  During 
this  stirring  process  clouds  of  air  bubbles  were  released  from  within  the  sediment,  and 
after  that,  the  sediment  texture  changed  completely,  becoming  much  "silkier"  to  the 
touch. 

The  reflected  signals  then  became  consistent,  as  can  be  seen  in  Figure  4.4,  where  pulses 
scattered  from  this  sediment,  after  each  stirring,  are  displayed.  The  inhomogeneities 
found  by  Chang  in  his  experiment  appear  to  have  been  caused  by  occluded  gas  within 
the  sediment.  Thereafter,  each  sample  of  sediment  being  studied  was  stirred 
periodically,  to  remove  any  air  bubbles. 

2.  Measurement  Considerations 

The  measurements  made  on  the  aquarium  sand  with  this  new  tecnique  still 
revealed  some  variation,  mostly  in  the  last  part  of  the  pulse  tails,  as  can  be  seen  in  the 
group  of  five  waveforms  in  Figure  4.5  and  in  the  measured  decay  slopes  presented  in 
Chapter  5. 

Measurements  for  the  aggregate  were  repeated  to  determine  if  the 
inconsistencies  in  our  results  were  also  due  to  the  presence  of  gas  bubbles  '\\'ithin  the 
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Figure  4.4     Waveforms  of  pulses  scattered  from  fine  sand, 
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Figure  4.5     Waveforms  of  pulses  scattered  from  aquarium  sand. 
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Figure  4.6    Waveforms  of  pulses  scattered  from  aggregate. 
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sample  or  due  to  the  visible  physical  inhomogeneities  present.  This  time,  the  variations 
remained  (Figure  4.6)  although  when  the  sediment  was  deeply  stirred  only  a  few 
bubbles  were  observed  during  the  stirring.  A  conclusion  can  be  drawn  from  these 
experiments  with  the  different  sediments  :  The  variations  in  the  tail  shape  of  the  pulse 
increases  as  the  grain  size  of  the  sediment  increases. 

The  consistent  dependence  of  the  results  on  the  grain  size,  indicated  that  we 
may  have  been  working  with  grain  dimensions  too  large  relative  to  the  wavelength.  In 
fact  the  average  dimension  of  the  aggregate  is  5.3  mm  and  the  wavelength  is  8  mm, 
values  that  are  very  close.  One  of  the  assumptions  necessary  to  determine  the  volume 
reverberation  is  that  the  surface  scattering  is  negligible,  because  the  sediment  surface  is 
levelled  and  smoothed.  The  levelling  of  the  bottom  surface  was  made  by  using  a  device 
that,  rolling  on  the  tops  of  the  tank  walls,  produced  a  horizontally  levelled  bottom 
surface.  Even  with  levelling,  when  the  grain  sizes  are  relatively  big  there  is  always  a 
gap  between  adjacent  grains  that  destroy  the  smoothness  and  produces  the  effect  of  a 
rough  surface. 

A  simple  experiment  was  performed  to  determine  the  importance  of  surface 
scattering.  A  Styrofoam  plate,  to  simulate  a  pressure  release  interface  was  weighted 
down  on  the  sediment  surface  and  covered  by  a  thin  layer  of  sediment,  and  the 
scattered  pulses  were  analysed.  In  this  case,  any  changes  in  the  pulse  waveforms  would 
result,  not  from  volume  reverberation,  but  from  surface  scattering. 
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Figure  4.7    Waveform  of  pulse  reHccted  from  a  thin  layer  of  aggregate. 

It  was  sufficient  to  observe  the  waveform  on  the  oscilloscope  to  fmd  that  the 
decay  changed  drastically  each  lime  the  transducer  was  slightly  moved.  Merely  movmg 
some  grains  from  the  thin  layer  caused  the  received  pulse  to  change  completely.  The 
piilbcs  (Figure  4.7  to  Figure  4.9)  showed  many  of  the  same  features  as  were  observed 
during  the  volume  reverberation  measurements.  It  was  concluded  that  surface 
scattering  was  a  significant  factor  in  reflection  from  sediment  when  the  grain  size  was 
comparable  to  the  wavelength  of  the  sound. 

Performing  the  same  experiment  with  Aquarium  sand  led  to  a  similar 
conclusion  that  surface  scattering  was  playing  an  important  role.  Although  the  pulse 
uaveforms  did  not  show  large  variations,  the  values  for  the  decay  slopes  showed  a 
significant  spread,  indicating  the  presence  of  surface  scattering. 

To  analyse  the  extent  of  the  inlluence  of  inhomogeneities  on  the  results 
obtained  in  the  measurements  performed  with  the  three  previously  indicated  sediments, 
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Figure  4.8     Waveform  of  pulse  reflected  from  a  thin  layer  of  aggregate. 

another  set  often  observational  nins  was  made  using  glass  beads.  It  was  expected  that 
the  reflected  pulses  could  be  more  consistent  due  to  their  homogeneity.  Ihis  did  not 
happen,  however,  and  Figure  4.1 1  shows  five  dilTerent  waveforms  reflected  by  the  glass 
beads.  They  present  approximately  the  same  features  as  the  pulses  reflected  by  the 
aquarium  sand  whose  average  grain  size  is  close  to  that  of  the  beads  dimensions, 
showing  once  more  the  influence  of  the  grain  size  in  its  scattering  properties. 
3.  Data  Processing 

The  computer  programs  used  in  this  project  were  designed  and  tested  by 
Chang  (1986).  They  were  used  to  obtain  the  average  slope  for  each  sample: 
THFSIS3  is  a  program  designed  to  convert  and  normalize  the  waveform  data  obtained 
from  the  Nicholet  digital  oscilloscope.  Since  the  pulses  contain  positive  and  negative 
values,  these  are  squared  and  a  series  of  fifty  waveform  samples  are  averaged  to  give  an 
accurate  representation  of  the  waveform  envelope.    Because  of  jitter  in  the  trigger 


40 


Figure  4.9    Waveform  of  pulse  reflected  from  a  thin  layer  of  aggregate. 

signal,  each  "bin"  of  the  scope  will  sample  the  waveform  at  various  times  during  a 
cycle,  and  the  resulting  curve  will  be  a  representation  of  the  envelope  of  the  pulse. 

The  result  can  be  seen  in  Figures  4.12  to  4.14  where  four  tails  of  the  pulses 
reflecicd  from  the  three  different  sediments  are  plotted. 

M.AXIMLM  is  a  program  used  to  choose  the  local  maxima  for  the  curses  obtained 
from  Til  ESI  S3  stored  on  disk.  The  square  root  of  these  maxima  is  then  taken  and 
stored  in  the  same  disk.  The  result  is  shown  in  Figures  4.15  to  4.17  whn.h  show 
respectively  the  envelopes  of  the  pulse  tails  shown  in  Figure  4.12  to  Figure  4.14. 
PLOT3  is  used  to  produce  the  plots  of  the  waveform  envelopes  and  to  apply  a  least 
square  method  to  fmd  the  best  fitted  straight  line,  its  slope  and  correlation  coelTicient. 
AVER.'XGF.  The  measurement  of  the  slope  of  the  tail  was  very  sensitive  to  any 
change  in  the  sediment  sample,  and  to  get  a  consistent  result  ten  sets  of  measurements 
were  performed,  and  this  program  is  applied  t.o  compute  a  chain  average  in  the 
following  way: 
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Figure  4.1 1     Waveform  of  pulses  scattered  by  glass  beads. 
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Figure  4.12    Tail  waveform  of  pulse  reflected  from  fme  sand. 
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Figure  4. 13    Tail  waveform  of  pulse  reflected  from  aquarium. 
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Figure  4.14    Tail  waveform  of  pulse  reflected  from  aggregate. 
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Figure  4.15    Tail  envelope  of  a  pulse  reflected  from  line  sand. 
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Figure  4.16    Tail  envelope  of  a  pulse  reflected  from  aquarium. 
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Figure  4.17    Tail  envelope  of  a  pulse  reflected  from  aggregate. 
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Each  averaging  step  is  divided  by  the  number  of  terms  envolved  so  that  each 
data  set  has  the  same  weight.  The  DAT  10  is  then  the  file  that  contains  the  averaged 
maxima  that  plotted  in  Figure  4.18  shows  a  much  more  regular  envelope  than  the 
envelopes  of  the  local  maxima  plotted,  for  example,  in  Figure  4.12. 


2 

^ 

CORRELATION  COEFF- 

9123 

1 

- 

SLOPE   -      -.0275 

0 

^ 

. 

::5i^ 

a" 

-1 

^^^ 

->^^ 

UJ 

< 
»— 
_J 
a 

> 

V 

z 
_J 

-2 
-3 
-4 

-5c 
c 

u 

1 

1 

1 

TIME    (MICRO-SEC) 
1         .         1         .         1         . 

1 

1 

3                         O 

3                         C\J 

n                  IT) 

ID 

O 

in 

o                 o                a 

CD                       O                       f\) 

m                 (O                ID 

o 

ID 

O 
(D 

Figure  4.18     Final  envelope  of  an  averaged  pulse. 
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V.  RESULTS 

The  decay  constants  of  the  trailing  edge  of  sound  pulses  reflected  off  the  bottom 
may  be  used  to  determine  the  volume  reverberation  of  bottom  material.  However,  the 
volume  reverberation  obtained  is  quite  sensitive  to  inhomogeneities  in  both  small-scale 
bottom  topography  and  sediment  composition. 

In  fact,  even  very  homogeneous  surfaces,  such  as  the  air  -water  interface  led  to 
slope  small  fluctuations  between  successive  runs.  For  this  reason,  averages  of  several 
runs  were  used  to  determine  the  volume  reverberation  level.  The  use  of  an  average  of 
several  runs  becomes  even  more  important  when  the  sediment  grain  size  increases. 

In  the  results  presented  next,  this  fact  is  very  evident.  In  Figure  5.1,  the  decay 
constants  obtained  for  the  different  sediments  show  big  fluctuations  mainly  for  the 
aggregate  and  aquarium  sand. 

A  proportionality  between  the  sediment  grain  size  and  the  amount  of  volume 
reverberation  was  expected  as  suggested  by  Clark,  Proni,  Seem  and  Tsai  (1984).  This 
was  observed  by  the  averaged  fmal  values  of  the  slopes  obtained  for  the  four  different 
sediments. 

A.       FINE  SAND 

For  the  fine  sand,  eight  measurements  were  made  and  the  individual  slopes  are 
listed  in  Table  2.    For  each  slope,  the  correlation  coefficient  was  also  computed. 

Despite  the  relative  stability  of  the  individual  slopes,  the  chain  average  was 
applied  to  the  8  measurements  and  the  results  are  listed  in  Table  3. 

The  running  average  of  all  eight  data  sets  is  ^  =  3.51  x  10  Np/^is  and  when 
compared  to  the  decay  constant  of  the  purely  reflected  (water-air  interface)  pulse,  a  = 
3.77  x  10"^  Np/jis,  shows  a  small  amount  of  volume  reverberation.  This  is  justified  due 
to  the  small  dimensions  of  the  grain  size,  (0.3  mm)  compared  with  the  8  mm  of  the 
wavelength,  making  a  a/X  =  0.375  which  can  be  considered  inside  the  Rayleigh  regime. 
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TABLE  2 

FINE  SAND  INDIVIDUAL  SLOPES 

DATA 

SLOPE  X  10^ 

CORRELATION 

SET 

(Np/Ms) 

COEFFICIENT 

DATA41 

-3.91 

-0. 9731 

DATA42 

-3.  42 

-0. 9162 

DATA43 

-3.  48 

-0. 9611 

DATA44 

-3.  58 

-0. 9682 

DATA45 

-3.  66 

-0.9575 

DATA46 

-3.26 

-0. 9526 

DATA47 

-3.  30 

-0. 9485 

DATA48 

-3.53 

-0.9871 

TABLE  3 

FINE  SAND  AVERAGE  SLOPES 

NO.  SETS 

SLOPE  X  10^ 

CORRELATION 

AVERAGED 

(Np/ns) 

COEFFICIENT 

2 

-3.  70 

-0. 9619 

3 

-3.66 

-0. 9715 

4 

-3.  61 

-0. 9721 

5 

-3.62 

-0. 9765 

6 

-3.62 

-0. 9703 

7 

-3.  50 

-0.9727 

8 

-3.51 

-0. 9757 
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B.       AQUARIUM  SAND 

A  series  of  ten  sets  of  measurements  were  run  for  this  sediment  (Table  4).  The 
individual  slopes  present  a  big  variation  (Figure  5.1).  Chain  averaging  gave  the  results 
in  Table5. 


TABLE  4 

AQUARIUM  INDIVIDUAL  SLOPES 

DATA 

SLOPE  X  10^ 

CORRELATION 

SET 

(Np/Ms) 

COEFFICIENT 

DATA51 

-3.00 

-0. 7717 

DATA52 

-2.  78 

-0.  9424 

DATA53 

-2.  68 

-0.7402 

DATA54 

-1.  88 

-0. 6269 

DATA55 

-2.43 

-0. 9242 

DATA56 

-3.42 

-0. 9557 

DATA57 

-3.03 

-0. 9777 

DATA58 

-2.35 

-0. 8495 

DATA59 

-3.  34 

-0. 8564 

DATA60 

1 

-3.  77 

-0. 9877 

The  final  decay  constant  Han  ~  2. 75  x  10""'  Np/)is  is  now  quite  different  from 
the  reference  value  of  the  water-air  interface.  This  difference  results  from  the  effect  of 
volume  reverberation  in  the  sediment  and  also,  as  shown  in  a  later  measurement,  from 
surface  scattering  due  to  the  relatively  big  dimensions  of  its  grains  compared  with  the 
wavelength  (  a/X  =  0.25  ). 
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TABLE  5 

AQUARIUM  AVERAGED  SLOPES 

NO.  SETS 

SLOPE  X  10^ 

CORRELATION 

AVERAGED 

(Np/Ms) 

COEFFICIENT 

2 

-3.05 

-0. 9634 

3 

-2.93 

-0. 9506 

4 

-2.53 

-0.8540 

5 

-2.  42 

-0.8782 

6 

-2.  67 

-0. 9003 

7 

-2.68 

-0. 9157 

8 

-2.  68 

-0. 8982 

9 

-2.  72 

-0. 9046 

10 

-2.  75 

-0. 9123 

C.  AGGREGATE 

Ten  measurements  were  performed  using  the  aggregate  with  the  resulting 
individual  slopes  hsted  in  Table  6.  The  variation,  in  the  presented  individual  slopes  is 
considerably  larger  than  the  variations  for  the  aquarium  sand  (Figure  5.1),  which  is  to 
be  expected,  since  the  aggregate  shows  many  inhomogeneities.  To  make  these  results 
more  understandable  the  chain  average  was  run  again  (Table  7). 

The  final  decay  constant  r|^g  =  1.84  x  lO"-^  Np/^s,  when  compared  v\ith  those 
for  the  other  sediments,  shows  an  even  larger  amount  of  reverberation  in  accordance 
with  the  proportionality  between  reverberation  level  and  size  of  individual  scatterers. 

D.  GLASS  BEADS 

Although  the  glass  beads  appear  homogeneous  in  size,  shape,  and  material,  the 
individual  slopes  show  considerable  variations  (Figure  5.1)  with  the  numerical  results 
listed  in  Table  8.   The  corresponding  averaged  values  are  presented  in  Table  9. 
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TABLE  6 

AGGREGATE  INDIVIDUAL  SLOPES 

DATA 

SLOPE  X  10^ 

1 
CORRELATION 

SET 

(Np/ns) 

COEFFICIENT 

DATA 51 

-2.  10 

-0. 9764 

DATA62 

-3.  54 

-0.9385 

DATA63 

-2.  40 

-0.8392 

DATA64 

-2.  71 

-0. 9772 

DATA65 

-1.  36 

-0.9322 

DATA66 

-2.  03 

-0. 8901 

DATA67 

-2.  09 

-0. 8187 

DATA68 

-1.55 

-0.9313 

DATA69 

-1.  13 

-0. 6211 

DATA70 

-1.84 

-0.9512 

The  final  decay  constant  for  the  glass  beads  is  t^-  =  3.06  x  10'^  Np/^s  which  is 
consistent  with  the  final  slopes  of  the  other  sediments,  considering  their  relative  grain 
sizes. 

E.       OVERALL  RESULT 

Relative  values  of  reverberation  level  for  the  different  materials  are  shown  in 
Figure  5.1  The  averaged  values  of  the  slopes  are  marked  along  the  y-axis  for  the 
different  materials.  Plotting  them  against  grain  size,  the  corresponding  best  fitted 
straight  line  is  defined  by: 

n  =    3.62  -  0.0035a 

The  value  of  the  average  grain  size  of  the  sediment  can  then  be  estimated  by  : 

a  =   10.34  -  n;0.0035 
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TABLE  7 

AGGREGATE  AVERAGE  VALUES 

NO.  SETS 

SLOPE  X  10^ 

CORRELATION 

AVERAGED 

(Np/fis) 

COEFFICIENT 

2 

-2.66 

-0.9476 

3 

-2.59 

-0. 9474 

4 

-2.  62 

-0. 9595 

5 

-2.  05 

-0. 9890 

6 

-2.05 

-0.9736 

7 

-2.  05 

-0. 9730 

8 

-2.00 

-0. 9756 

9 

-1.  85 

-0. 9653 

10 

-1.  84 

-0. 9719 

where  s  is  in  Np/^s  and  a  is  in  mm.  This  relationship  gives  only  a  rough  estimate  of 
the  sediment  grain  size  since,  considering  the  standard  deviations  in  the  slopes  of  the 
extreme  data  points,  the  value  for  a  can  be  estimated  within  ±  2.8  mm.  for  a  slope  of 
2.5  X  10  Np'^s.  Also  this  relation  was  established  based  on  measurements 
performed  using  a  transducer  with  damping  corresponding  to  a  decay  constant  a  = 
3.77  X  10'^  Np,^s  at  180  kHz.  Other  transducers  will  have  different  decay  constants 
and  further  measurements  should  be  performed  using  different  transducers  to 
investisate  the  influence  of  a  on  this  result. 
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TABLE  8 

GLASS  BEADS  INDIVIDUAL  SLOPES 

DATA 

SLOPE  X  10^ 

CORRELATION 

SET 

(Np/Ms) 

COEFFICIENT 

DATA71 

-3.48 

-0. 9896 

DATA? 2 

-2.92 

-0. 915? 

DATA? 3 

-2.99 

-0.9565 

DATA? 4 

-3.03 

-0.8251 

DATA? 5 

-2.  10 

-0.8131 

DATA? 6 

-3.22 

-0. 9286 

DATA?? 

-3.81 

-0.991? 

DATA? 8 

-3.40 

-0.9860 

DATA? 9 

-3.94 

-0.8614 

DATA80 

-2.92 

-0.9680 
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TABLE  9 

GLASS  BEADS  AVERAGED  SLOPES 

NO.  SETS 

SLOPE  X  10^ 

CORRELATION 

AVERAGED 

(Np/^s) 

COEFFICIENT 

2 

-3.  33 

-0. 9794 

3 

-3.  19 

-0. 9750 

4 

-3.28 

-0. 9646 

5 

-2.  86 

-0. 9421 

6 

-2.  84 

-0. 9693 

7 

-2.83 

-0. 9780 

8 

-2.94 

-0. 9829 

9 

-3.  01 

-0. 9789 

10 

-3.  06 

-0. 9805 
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Figure  5.1     Spread  of  individual  slopes. 
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VI.  CONCLUSIONS 

The  purpose  of  this  experiment  is  to  establish  a  relationship  between  acoustic 
volume  scattering  coefficients  and  sediment  grain  sizes.  In  the  experimental  design 
only  volume  reverberation  was  to  be  measured,  and  surface  scattering  was  to  be 
considered  negligible.  Surface  scattering  was  not  negligible  for  most  of  the  materials 
studied,  and  it  was  not  possible  to  achieve  the  initial  purpose. 

Nevertheless  some  conclusions  can  be  drawn: 

1.  It  was  found  that  sediment  grain  size  is  related  to  the  decay  constant  of  the 
tail  of  the  scattered  pulse  by  the  equation 

a  =  10.34  -  Ti/0.35  (6.1) 

This  equation  has,  however,  some  limitations  since  it  is  applied  to  the  decay  constant 
of  the  specific  transducer  used  and  to  the  frequency  of  180  kHz.  It  can  be  used, 
however,  if  the  transducer  is  calibrated  to  a  decay  constant  a  =  -3.77  x  10" •^  Np'^s. 

2.  For  wavelengths  of  the  same  order  of  magnitude  as  the  average  grain  size, 
surface  scattering  is  an  important  component  of  the  reverberation  level,  even  if  the 
sediment  surface  is  perfectly  levelled  and  smoothed,  and  depends  also  on  the  sediment's 
grain  size. 

3.  The  frequency  used  and  materials  selected  for  the  experiment  correspond  to  a 
wide  range  of  ka  values,  from  0.12  for  iht  fine  sand  to  4.16  for  the  aggregate.  The 
decaying  slope  monotonically  decreases  with  increasing  grain  size  and  does  not  show 
any  strange  behavior  in  the  region  of  ka  =  1. 

4.  Inhomogeneities  in  sediment  composition  are  not  responsible  for  the  variation 
in  the  decaying  slopes  obtained,  at  least  not  totally,  since  measurements  made  on 
homogeneous  material  like  glass  beads  show  the  same  spread  in  results  as  the  aquarium 
sand  of  similar  size. 

Further  study  of  the  volume  reverberation  from  ocean  sediments  is 
recommended,  as  is  the  use  of  different  combinations  of  frequency  and  grain  sizes, 
keeping  in  mind  that  wavelengths  of  the  same  order  of  magnitude  as  the  grain  size 
should  not  be  used  so  that  surface  scattering  may  be  neglected. 
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A  theoretical  search  must  be  done  to  investigate  the  influence  of  transducer 
characteristics  on  the  decay  constant,  to  venfy  the  extent  of  the  validity  of  Equation 


6.1 


59 


APPENDIX 
DATA  COLLECTING  AND  PROCESSING  SEQUENCE 

In  this  appendix  a  complete  sequence  of  the  data  collecting  and  processing  for 
the  fine  sand  is  presented.  Eight  data  sets  (DAYAU)  of  50  pulses  each  were  collected 
for  this  sediment  and  averaged  using  the  program  THF:SIS3.  These  data  sets  (files 
DATA41  to  48)  are  graphically  represented  in  Figures  1  to  8  where  the  squared 
voltages  are  plotted  against  time  (program  PLOT). 


Figure  1     Pulse  waveform  of  DATA41. 

The  decaying  slopes  for  these  pulses  are  then  determined  using  the  program 
MAXIMUM  that  selects  the  points  of  local  maximum  voltages.  A  list  of  these  points 
on  the  pulse  tail  for  the  first  data  set  (DArA4l)is  presented  next. 
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LISTING  OF  LOCAL  MAXIMA  FOR  DATA41 
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The  corresponding  waveform  envelope  is  plotted  in  Figure  9.  It  is  obtained  by 
connecting  the  points  of  maximum  voltage.  The  best  fitted  straight  line  and 
corresponding  slope  are  then  determined  (program  PLOl  3). 

The  individual  slopes  show  big  spread  (see  Chapter  five),  and  they  were  applied 
in  the  averaging  process.   The  eight  data  sets  were  averaged  as  follows: 
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Figure  2     Pulse  waveform  of  DATA42. 
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Figure  3     Pulse  waveform  of  DATA43. 
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Figure  5     Pulse  waveform  of  DATA45. 
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Figure  6     Pulse  waveform  orDAlA46. 
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Figure  7    Pulse  waveform  of  DATA47. 
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Figure  S     Pulse  waveform  of  DA  rA48. 


D/\T42  is  the  average  of  DATA41  and  DATA42.  From  DAT42  to  DAT48  the 
variations  presented  in  the  individual  waveforms  are  filtered  because  the  nuniber  of 
terms  involved  in  the  average  increases.  DAT48,  the  result  of  the  averaging,  may  be 
considered  the  representative  waveform  for  the  fine  sand. 

Next,  programs  .MAXIMUM  and  PL0T3  are  applied  to  the  averaged  DAT#  sets 
which  corresponding  plots  are  presented  in  Figures  10  to  17. 

The  envelope  from  DATA41  to  DAT48  gets  smoother  and  the  value  of  the  decaying 
slope  stabilizes  along  the  averaging  process  although  the  differences  are  small  since  the 
individual  waveforms  are  not  very  different  for  this  sedhnent. 
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